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Correlations are important tools in the characterization of quantum fields, as they can be used to 
describe statistical properties of the fields, such as bunching and anti-bunching, as well as to perform 
field state tomography. Here we analyse experiments by Bozyigit et al. [T] where correlation functions 
can be observed using the measurement records of linear detectors (i.e. quadrature measurements), 
instead of relying on intensity or number detectors. We also describe how large amplitude noise 
introduced by these detectors can be quantified and subtracted from the data. This enables, in 
particular, the observation of first- and second-order coherence functions of microwave photon fields 
generated using circuit quantum-electrodynamics and propagating in superconducting transmission 
lines under the condition that noise is sufficiently low. 



I. INTRODUCTION 



Field correlations are widely used in the characteriza- 
tion of classical and quantum fields [21 [3] . A particular 
set of correlations used for such purposes are the coher- 
ence functions of a field, as described by Glauber [JHS]- 
These functions can be used to quantify the ability of 
a field to interfere with itself, as well as to demonstrate 
features of quantum fields which cannot be reproduced 
in a classical system. One of the most famous of these 
quantum phenomena is known as anti-bunching [71 [S], 
and it is frequently used to characterize single-photon 
sources in the optical regime ^-[ii . Over the recent 
years Josephson-junction based superconducting circuits, 
resonators and transmission lines have emerged as a plat- 
form for performing quantum optics experiments in the 
microwave regime [Tm24j . While in the optical regime 
coherence functions are usually measured using an inter- 
ferometer where photon number detectors are used, in 
the microwave regime, linear detectors {i.e. field quadra- 
ture measurements) are ubiquitous due to the difficulty 
of building reliable photon number detectors. This raises 
the question of how to measure field correlations us- 
ing linear detectors. This paper answers this question 
and describes the theory behind the recent experiments 
performed by Bozyigit et al. where correlations of 
a propagating microwave field are measured using only 
linear detectors, instead of intensity detectors. While 
the discussion here focuses on the measurement of first- 
and second-order coherence functions of microwave fields, 
the analysis can be applied to any correlation of field 
operators. We note that the measurement of correla- 
tion functions of propagating microwave fields using non- 
linear (i.e. square-law) detectors was theoretically stud- 
ied in Ref. [25], under the assumption of negligible cor- 
relation in the noise added by the detection chain. In 
practice, these correlations turn out to be important and 
are discussed here. Recent work by Menzel et al. [26] and 
Mariantoni et al. [27^ is in a similar direction to the work 
presented here. 



The paper is organized as follows. Section [TT] gives 
a brief review of coherence functions and how they are 
measured with non-linear detectors. Section lllll describes 
how field correlations, and in particular coherence func- 
tions, can be measured using linear detectors. Section IV 
describes the effects of noise in the experiments, and fi- 
nally Section [V] describes how the experimental setup can 
be simplified in circuit QED experiments. 



II. COHERENCE FUNCTIONS 

The meaning of coherence of a field in a single fre- 
quency mode, with corresponding annihilation operator 
a, can be understood by considering interference experi- 
ments which use the field leaking out of this mode. Using 
a double slit, the field can be made to travel two pathways 
of different lengths which terminate at a single point-like 
photon detector, as depicted in Fig. [ija) . The combined 
field that impinges on the detector is made up of fields 
originally emitted at times t and t + t (which depend 
on the lengths of the paths), so that the observed field 
intensity at the detector is the sum of the intensities of 
the two fields plus an interference term which depends 
on {a^t)a{t + t)) [2]. Interference effects can only be 
observed if this correlation is non-zero. It is therefore 
natural to define 



G^^\t,t + T) = {a\t)a{t + T)), 



(1) 



which is called the first-order coherence function [4], as a 
measure of the emitted field's potential to interfere with 
itself - in other words, a measure of the coherence of the 
field. One may also consider 



(2) 



for some time interval T in order to obtain an expression 
that depends only on the time difference between the 
two paths. If G'-^^ (r) = 0, no interference effects can be 
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FIG. 1: (Color online) Experimental setups illustrating dif- 
ferent degrees of optical coherence with intensity detectors. 
The red lines represent quantum fields, and the black lines 
represent measurement records. 



observed for a path difference of cr, where c is the speed 
of light. 

The measure of coherence that is most often used to 
distinguish classical fields from quantum fields is the 
second- order coherence function given by 

G(2) {t, t + T) = {d^ {t)a^ (t + T)d{t + r)a(t)) , (3) 

or by the integrated version 



(4) 



The canonical experiment which gives the physical in- 
terpretation of G'^-' is one with a single light source and 
two point-like detectors, such that the field takes a time t 
and < + r respectively to reach each detector, as depicted 
in Fig. [ijb). In that case the correlation between the 
detected intensities is given by G^^^(i, t + r). 

For classical fields, where the field operator in the ex- 
pressions above are replaced by c-numbers, one finds that 
|G(^'(0)| > \G^'^\t)\, while there are quantum states of 
the field that yield |G(2)(0)| < |G(2)(t)| for r 9^ 0, a phe- 
nomenon known as anti-hunching [318]. The canonical 
examples of anti-bunched field states are single photon 
states and squeezed states. In the case of pulsed exper- 
iments - where the light field state is prepared with a 
repetition period of ip - one writes instead that classical 
fields obey |G(2)(0)| > \G^'^\ktp)\, and that some quan- 
tum states of the field yield |G(2)(0)| < \G'^^\ktJ)\ for 
/c 7^ 0. Only pulsed experiments will be considered in the 
remainder of this paper, the generalization to continuous 
experiments being straightforward. 



A. Standard experimental setups 

As illustrated in Fig. [2] we consider experiments where 
the source is a single mode of a cavity coupled to a trans- 
mission line via a leaky mirror, a situation typical of cav- 
ity QED [251 - 150] . In circuit QED for example, arbitrary 
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FIG. 2: (Color online) Standard experiments for the observa- 
tion of G<^' and G'^' using intensity detectors: (a) a Mach- 
Zender interferometer with a variable delay r and a variable 
phase shift ip, and (b) a Hanbury Brown and Twiss (HBT) 
interferometer with a variable delay r. The light-blue com- 
ponents are balanced beam-splitters. The cavity is taken to 
be one-sided, with one mirror being perfectly reflective. 



superpositions of a single photon and vacuum can be pre- 
pared in the dispersive regime via Purcell decay |17j or by 
strong coupling to a qubit brought into resonance with 
the cavity J^, although details of the state preparation 
are not important for the remainder of the discussion. 
The harmonic field in the cavity is associated with an 
annihilation operators d with the usual same-time com- 
mutation relation [d,a)] = 1. Using input-output the- 
ory [3l f3Tti33] . one can show that d is related to the modes 
of the transmission line via 



bout{t) = y/i^d{t) - bin{t), 



(5) 

where Kb is the rate at which photons leak out of a, and 
the input and output fields are given by 
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da;e"*'^(*-*i)6(ii,a;) (7) 



for transmission line modes b(t,Lu) at times to < t < ti, 
and correspond to fields propagating towards or away 
from the cavity. The commutation relations of the input 
and output fields are given by 



[bUt),bl{t + T)] = [6,,tW,&Lt(i + ^)] = '5(t). (8) 

These definitions lead to an equation of motion for d in 
the interaction frame to be given, for a one-sided cavity, 
by 



d{t) 



-^fait) 



Kbhnit). 



(9) 
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From Eq. ^ it is clear that the correlations of 5out are 
proportional to the correlations of a when &g is prepared 
in the vacuum state. The remainder of the discussion 
will focus on the observation of the coherence functions 
of the output field &out only, as they can be taken to be 
equivalent to the correlation functions of a. The "out" 
subscript will also be dropped when it is clear form the 
context. 

The state of the cavity field is taken to be prepared at 
times t = ktp for integer k and repetition period tp, and 
allowed to decay via the leaky mirror as described above. 
The repetition period is chosen to obey tp ;» 2n/Kh so 
that the cavity can be taken to be in equilibrium at the 
time of the next preparation of the cavity field. 

When working with photons in the optical frequencies, 
G'-^-' is usually observed using a Mach-Zender interferom- 
eter with a variable delay of r in one of the branches [31] , 
as depicted in Fig. [2ja). The difference between the in- 
tensities in the photo-current detectors can yield the real 
or the imaginary part of G^^\ depending on the phase 
shift if in the lower branch. The standard approach to 
the observation of the G'^-* is to use a Hanbury Brown 
and Twiss (HBT) interferometer [34| , which is illustrated 
in Fig. [2|b). In order to observe G^^\ one simply mea- 
sures the correlations between the photo-currents of the 
two detectors. 

Both these setups rely on field intensity detectors, 
which give information about the number of photons, 
and thus can be modeled by non-linear quantum opti- 
cal interactions^. Low- noise intensity detectors for opti- 
cal fields are common, and although non-linear detectors 
have been demonstrated in the microwave regime |17l I35j 
(albeit with higher noise levels than in optics), the main 
motivation for this paper is to illustrate how the coher- 
ence functions of microwave fields in circuit QED may be 
measured through the use of linear detectors only. 



III. LINEAR DETECTORS 

Field quadrature measurements of microwave signals 
is a standard technique which has been applied 
very successfully to quantum electrical circuits in the 
recent years to demonstrate, for example, new regimes 
of cavity QED [T3], high-contrast detection of qubit 
states |37], photon states [T5], and nanomechanical os- 
cillator states [38^. Since field quadrature operators are 
fundamentally different from number operators, different 
experimental setups are required in order to measure the 
coherence functions G'^' and G^'^\ Grosse et al. [SHj have 
demonstrated how a HBT interferometers can be mod- 
ified to measure using field quadratures instead of 



Linearity in this sense refers to tlie representability of tlie Heisen- 
berg picture evolution by a linear transformation of creation and 
annihilation operators for all times. 
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FIG. 3: (Color online) Optical analog of an IQ-mixer as an 
8-port homodyne detector. The field f is fed into a balanced 
beam splitter along with vacuum Vr. The X quadrature of 
one of the beam splitter outputs is measured while the P 
quadrature of the other output is measured. The classical 
outcomes are the real and imaginary parts of the complex 
envelope Srit). 



intensity measurements. Here we analyse similar exper- 
iments [II, and consider generalizations and simplifica- 
tions which exploit features of circuit QED, while at the 
same time considering the large added noise due to the 
HEMT amplifiers currently required for measurement in 
this system. 

The details of the implementation of quadrature op- 
erator measurements in the microwave regime are differ- 
ent from the standard optical implementation. In par- 
ticular, homodyne detection in the microwave regime is 
performed via mixing instead of beam splitting |36| . For 
simplicity, we will however consider the optical analogues 
of the devices we discuss. Common non-idealities in the 
microwave regime, such as weak thermal states instead of 
vacuum inputs, can be treated straightforwardly by con- 
sidering different input states, and thus do not change 
the analysis significantly. 

The measurement of both quadratures of a propagat- 
ing field, realized in optics through 8-port homodyne |3D] 
or heterodyne detection, is performed by an IQ mixer in 
the microwave regime |36j . The symbol for the IQ mixer, 
and its description in terms of its optical analogue are 
depicted in Fig. [3j The input is any propagating quan- 
tum field with annihilation operator f, which may stand 
for any propagating field considered in this paper. The 
outputs are quadrature measurements of the superposi- 
tions of the f field with a mode Vj. in the vacuum state, 
where [f,{ij] = 0. These outputs are labeled Xi and P2 
to emphasize that the measurements are made on differ- 
ent commuting modes, and correspond to the in-phase 
component and the quadrature component of the mea- 
surement respectively. 

Finally, it is important to note that, for most circuit 
QED experiments, only averages of these quadratures 
over many realizations of the experiment are measured. 
Here, however, we are interested in experiments where 
the full time records of these quadratures are recorded, 
for each realizations of the experiments [I]. Based on 
these full records, any averages or correlation functions 
can be reconstructed, as is discussed in the next sections. 
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A. Complex envelope 

Given the two classical outputs Xi{t) and P2{t), it is 
useful to define the complex envelope Sr{t) of f as 



Sr{t)=X,{t)+iP2{t), 



(10) 



which is a random c-number due to the dependence on 
the measurement records of the quadratures. Noting that 

(Mt)) = ( ) - (Mt)) + (Mt)) (11) 



(P^it)) = -H ) = (Prit)) - (A,(0>, (12) 



one may write that {Sr{t)) — {Sr{t)) where the complex 
envelope operator Sr is defined by 



Sr{t) = f{t)+vl{t)^X,{t)+tP2{t). 



(13) 



In order to simplify the remainder of the calculations, it 
is convenient to define 5"^ in this manner instead of using 
the quadrature operators explicitly. 

Given that the mode Vr is in the vacuum state, the 
expression for the expectation values take simple forms. 
The presence of the vacuum mode Vr is indeed important 
as it leads to the commutation relation 



(14) 



implying that Sr is normal and therefore diagonalizable. 
Since Sr is described by the sum of the commuting oper- 
ators Xi and P2, its eigenvalues are given by the sum of 
the eigenvalues of these operators for any fixed eigenvec- 
tor. This corresponds to the measurement record Sr of 
Sr being simply the sum of the measurement records Xi 
and P2, as claimed earlier. Note that this does not imply 
that both quadratures of r can be measured simultane- 
ously without back-action. 

Since Sr and S^ commute at all times, arbitrary cor- 
relations of these operators, like the correlations of their 
measurement records, do not depend on operator order- 
ing. Therefore, 

((5;)"^;) = {(sir^s:) = {(f^+v^rif+vir), (is) 

independently of the ordering of the terms. However, in 
order to reduce these expressions to a correlation function 
of f alone, one must rewrite the expression such that the 
Vr modes are in normal ordering in order to immediately 
evaluate the expectation values, leading to 

((5;)™5;') = {{Sirs';) = (f"(ft)'"), (le) 

or in other words, correlations of the complex envelope of 
a field correspond to anti-normally ordered correlations of 
the field operator, under the assumption that the Vr mode 
is prepared in the vacuum. In this case the measurement 



-V 



IQ 



■Sc{t) 



IQ 
mix 



Sd{t) 

FIG. 4; (Color online) Hanbury Brown and Twiss interfer- 
ometer with complex envelope measurement used to measure 
the first-order coherence function G'^'. The output of the 
cavity field is separated by an on-chip beam splitter and the 
resulting fields measured by an IQ-mixer. 



of Sr is described by the Husimi-Kano Q function which 
is known to give access to anti-normally ordered same- 
time correlations [IDHH]. Similar results hold for multi- 
time correlations. 

It is important to note that, while with this approach 
arbitrary correlations can be evaluated, the number of 
statistical samples needed to obtain a desired precision 
in the estimate grows as the noise power raised to the 
desired correlation order (see Appendix |A] for details) . 
In practice, this limits the order of the correlations mea- 
sured with current amplifier noise levels due to the large 
number of repetitions of the experiment needed to obtain 
reasonable error bars. Use of quantum limited amplifiers 
would greatly improve the situation |2TJ [211 US] ■ 



B. G^^^ observation 

As depicted in Fig. [41 with IQ-mixers, the first order 
correlation function can be measured from the out- 
puts of a HBT interferometer. Because of the unitary of 
the beam splitter and the presence of a vacuum port, the 
complex envelope operators of the outputs labeled c and 
d commute. 

The auto-correlation of one of the complex envelopes, 
say Sc{t), is given by 

rW(i,t + r) = {Sl{t)S,{t + r)) = 6{r) + ^{bHt)b{t + r)), 

(17) 

while the cross-correlation between the complex en- 
velopes is 

rW(t,t+r) = {Sl{t)S,{t + r)) = ^{b^t)b{t + r)), (18) 

where we have used the fact that the expectation values 
of all the vacuum modes are zero. Thus the first-order 
coherence function G^^^ of the b field is immediately ac- 
cessible from cross-correlations of the complex envelopes 
in a modified HBT interferometer via 

GW(i,t + r) = 2rW(t,t + r)-25(r), 

.(1),, . . ^ (19) 



2r'^'> {t,t + T), 
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up to non-idealities, such as amplifier noise, which will 
be treated later. The exact expressions for G*^^' of the 
states prepared in Ref. [T] are given in Appendix [B| 

Although the divergence of the S functions may ap- 
pear problematic, in reality due to the finite bandwidth 
of the experiments these delta functions are replaced by 
smooth bounded functions, while the coherence functions 
are distorted by a convolution kernel which preserves the 
relative heights of the peaks in the experiment. This re- 
sults in the filtered correlation functions 



4l^M + /cff(T), 



(20) 



where Ggj''(r) = G(^^(T)*/off (r) and f^s is a function de- 
scribing the effective action of the filter (see Appendix [C| 
for details). 



C. G^^' observation 

The expressions needed to measure the second-order 
coherence function from the complex envelopes can be 



constructed by inspection from Eq. (13 1. Depending on 



which factors are taken to be complex conjugates or to 
be displaced in time by r, different correlations can be 
used to extract information about G^^\ One such choice 



IS 



+ hbHt)bHt + r)b{t + r)bit)), (21) 



so that G(2) can be obtained immediately via 

G(2) (t, t + T) = 4r(f ) {t, t + T)- 2G(i) (t, t)S{0) 
- 2G^^\t + T,t + t)5(0) - 4(52(0). 

Another choice that leads more directly to G'-^-' is 
Tf{t,t + T) = {Sl{t)Sl{t + T)Sa{t + T)SS)), 



so that 



^-{P{t)b'^{t + T)b{t + T)h{t)), 



G'^^\t,t + T)=ATf{t,t + T). 



(22) 



(23) 



(24) 



As described earlier, the divergence of the 5 functions is 
taken care of by filtering in a realistic experiment. The 
main distinction between these two approaches of mea- 
suring the second-order coherence functions is how they 
are affected by noise in the experiment, as is discussed in 
the next section. 



IV. REJECTION AND SUBTRACTION OF 
NOISE 



The amplitude of microwave signals in a supercon- 
ducting quantum circuit is small enough that ampli- 
fiers are essential for their observation, and so in a re- 
alistic experiment, the field is amplified before mixing. 
Using the Haus-Caves description of a quantum ampli- 
fier [331 1431 144| , an input operator c and an output op- 
erator Camp for a phase-preserving amplifier with gain gc 



are related by 



/5< 



(25) 



where he is an added noise mode. 

It is clear that if (7^ > 1 there will be added noise due 
to amplification, even at zero temperature. However, for 
thermal white Gaussian noise, one finds that all odd or- 
der moments vanish. As a result, the first moments of 
quadrature fields are not affected by this amplifier noise, 
just as they are not affected by vacuum noise. The contri- 
butions from other moments may be non-zero, however, 
and must be accounted for. For simplicity, we only con- 
sider the case of Gaussian white noise here, but similar 
results follow straightforwardly for general noise as long 
as the noise is independent of the inputs. Since the noise 
moments can be extracted from experimental data, the 
assumption of Gaussian noise is not essential. 

The noise modes from different amplifiers are taken to 
commute, but in general they may be correlated. While 
the noise is normally taken to come from the amplifi- 
cation [331 SSI IUj formally one may also take he to 
include thermal noise from other sources, such as the 
vacuum ports of the IQ-mixer and of the beam-splitter, 
with only minor modifications. Here he is taken to have 
a commutator [hc{t),h\{t -f r)] = 8{t) in order to pre- 
serve the bosonic commutation relations of the ampli- 
fied signals Camp, and auto-correlation {h\{t)hc{t + r)) = 
NcS{t). The noise sources are assumed to be indepen- 
dent of the inputs, so that [c, he] = [c, /ij] = 0, and 
(che) = (chl) = 0. The correlations between he and the 
noise mode hd from the other amplifier in the experiments 
described here is taken to be {he{t)h\{t + r)) = NedS{T) 
while {he{t)hd{t + t)) = 0. 

Using this noise model, one can calculate the different 
correlations F^^^^ using the amplified modes, resulting in 

riUmp(^>i + = jG^'Ht,t + r) + {Ne + 9c)5{r), 
r?Lp(^> t + r)= ^^G(i) (<, t + r)+ NedSir), 

(26) 

in the unfiltered case. 

Since the thermal noise in the amplifiers is indepen- 
dent of the inputs, a steady-state experiment with the 
input mode b in the vacuum state can be used to esti- 
mate the noise strengths and subtract the corresponding 
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terms from T^^^^ to obtain an estimate of G^^^ . When 
the noise cross-correlation Ned is expected to be zero 
or neghgible compared to the noise auto-correlations Nc 
and Nd, the approach to the estimation of G^-^^ based on 



^y^amp provides noise rejection without additional post- 
processing. 

The second-order coherence function for the amplified 
fields has similar properties. One finds 



r2'lmp(i,i + r) = ^-^G^^Ht,t + r) + |<5(0)[ff, + Nd]G'-'\t,t) + |5(0)[5c + Nc]G^'\t + r,t + r) 

+ ^^^N,d5{T)[G^^\t + T, t) + G(i)(t, t + r)] + [g,Nd + g^Qd + 9dNc]5\Q) + {hlit + T)h\{t)K{t)hd{t + r)), (27) 



while 



r;^'Lp(^. i + (i, t + r) + + T)h\{t)K{t)K{t + r)) 



Ned 



J (t)G(i) (t + T, t) + 5(0)G(i) {t + T,t + T)+ ^(0)G(i) (t, <) + (5(r)G(i) (i, t + r)! , (28) 
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FIG. 5: (Color online) The setup for the observation of co- 
herence functions using a two-sided cavity. 

where all odd moments of the noise modes where taken 
to be zero (if such an assumption cannot be made, simi- 
lar expressions involving the odd moments are easily de- 
rived but are omitted here for brevity). The recovery of 
the second-order coherence function from noisy signals is 
clearly more involved, but requires only the estimation of 
first-order coherence functions, as well as two and four- 
point noise correlations in an experiment where the input 
mode b is prepared in the vacuum. Since the filters are 

(21 

taken to be linear and time-invariant, Ggj is a scaled and 
distorted version of G(2), preserving the relative heights 
of the peaks, so that the non-classical properties of the 

(2) 

field can still be verified. Once again we see that ^^^^^p 
provides a more direct estimation of G^^^ by rejecting 
contributions from uncorrelated noise up to four-point 



i-SIDED CAVITIES 

Strictly speaking, the beam splitter is not necessary 
for the observation of the coherence functions described 
above. If one considers a two-sided cavity, illustrated 
in Fig. [Sj the correlations between the cavity outputs 
behave in a manner similar to the outputs of the beam 
splitter in the HBT interferometers. In particular, using 
causality as well as the boundary conditions of the input 
and output fields of the two-sided cavity. Appendix [P] 
shows that 

[V,t(i),cLt(i')]=0, (29) 

where Cout is defined in a manner analogous to 6out , with 
a mirror leakage rate Kc^ and an amplifier with gain 
being applied before mixing and measurement. It is thus 
possible to measure the complex envelopes of the two 
cavity outputs and calculate the correlations in the same 
manner as in the modified HBT setup without the need 
for an additional beam splitter. This can lead to sim- 
pler and smaller experimental setups, as beam splitters 
in the microwave regime can occupy a significant area in 
coplanar devices. 

Calculating the correlations using the two cavity out- 
puts 6out and Cout one finds 
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.(1) 



r)i;Lp(^' t + r) = (i, t + r) + iVbe(5(r), 

+ V%g~cV^cNic5{T) [G(i) (i + r, i) + G(i) (i, t + r)] 

+ [gbN, + gbgc + 9cNb]S^iO) + {hl{t + T)h\{t)}n,{t)h,{t + r)) 



[g, + N,]5{Q)gbKkG^^\t + T,t + T) 



.(2) 
/3,amp 



gbgcKbKcG'^^^it, t + T) + {hl{t + T)hl{t)hc{t)hc{t + r)> 



:7^^fcA^7Vb, {,5(r) [G(i) (i + r, i) + G^^) (i, t + r)] + 5(0) [G^^^ (< + r, i + r) + G^^^ {t, t)] } 



(30) 
(31) 

(32) 
(33) 



where G^^-* and G*-^-* are now the coherence functions of 
the cavity field a instead of the cavity output fields, lead- 
ing to the introduction of additional factors which depend 
on the cavity leakeage rates Hb,c- These expressions are 
directly analogous to Eqs. (p6|), (pTl, and (|28|). 



VI. SUMMARY 

We have analysed experiments for the measurement 
of field correlations using only field quadrature detec- 
tors and in the situation where the full record of many 
repetitions of the experiment are available. The combi- 
nation of the quadrature measurements into complex en- 
velopes gives direct access to anti-normally ordered field 
correlations. While re-ordering of the operators in the 
correlations and the use of phase-preserving amplifiers 
introduces additional noise into these measurements, we 
demonstrated that the noise can be accounted for and 
subtracted in order to reveal only the field correlations of 
interest. Although there are indications that the number 
of statistical samples scales exponentially with the order 
of the correlation function, the measurement of low order 
correlations is possible for current amplifier noise levels. 
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Appendix A: Statistical error on correlation 
function estimates 

In order to estimate the minimal number of repetitions 
of the experiment which must be performed to extract a 
given correlation function, consider the product of uncor- 
related Gaussian random variables with zero mean and 
identical variances cr^. These random variables corre- 
spond to the measurements of different outputs at steady- 



state after the cavity state has decayed, and the variances 
are given by the noise power of the measurement record 
(including vacuum noise). In order to illustrate the argu- 
ment, we consider real valued random variables Vi first, 
and generalize to complex valued random variables G^. 
Since these random variables are uncorrelated, it follows 
that {V1V2 ■ • ■ Vm) = 0. However, given a finite number 
of statistical samples, the sample average V1V2 • • • K« will 
deviate from zero due to statistical fluctuations. Signal 
features which are comparable with the typical size of 
these fluctuations cannot be reliably observed. As the 
typical size of these fluctuations decreases with the in- 
creasing number of repetitions, this is in principle not a 
fundamental problem. 

In order to estimate the number of samples needed for 
the reliable estimation of two-point correlations, consider 
the product of two Gaussian random variables. The char- 
acteristic function of this product is given by 



0([/) 



dvi dv2 du pi{vi)p2{v2)S{viV2 - u)e 

(Al) 
(A2) 



1 



The characteristic function of the average of R samples 
is given by 



(A3) 



Given some error e > and a number of repetitions 
R, the probability that the sample average V1V2 obeys 
— <^/2 < V1V2 < V2 is given by the integral of the inverse 
Fourier transform of (t>ii{U) over this range and simplifies 
to 



Ft {\V,V2\ < i2) - 



1 

2^ 



+ 00 



dUMU) 



sm 



(A4) 



which can be evaluated by numerically. Thus it is 
straightforward to calculate the number of repetitions R 
required to observe a feature larger than e with confi- 
dence Pr(|VVl^| < V2). 
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Another approach that provides a looser bound, but is 
more readily generalized to higher order correlations, is 
based on Chebyshev's inequality [15] . The variance of the 
product of independent random variables with zero mean 
is the product of the variances of each of the random 
variables. In the case of R samples of the product of m 
independent random variables Vi one finds that 



PT( ViV2---VrJ<i2) > 1-4 



(A5) 



Note that in order to obtain this bound no assumption 
was made about the form of distribution of the random 
variables, other than the fact that the random variables 
are independent. Solving for R one obtains the worst- 
case upper bound 

R < ,, (A6) 

which makes clear the exponential relationship between 
the order of the correlation and the number of samples 
needed to have a statistical error of less than «/2 with 
some fixed probability. 

In order to generalizing this to complex- valued random 
variables Ci - where the real and imaginary parts of Ci 
are independent with variance a, and the Ci are mutu- 
ally independent - simply consider the real and imagi- 
nary parts of the correlations separately. In that case, 
because a larger number of terms contribute to the real 
and imaginary parts of the correlation, the variance has 
a larger bound, and one finds 



R< 



g2 Pr (error) ' 



(A7) 



where Pr(error) is the probability that the absolute value 
of the real or imaginary parts of C1C2 ■ ■ • Cm are greater 
than '=/2. 

There is no indication that taking into account the 
Gaussian statistics of the random variables leads to bet- 
ter scalings. Thus the ratio of the number of statistical 
samples needed to estimate G^^-* vs. G'^^' for some fixed 
noise variance and desired accuracy is at worse propor- 
tional to the noise power in the experiments. As a result 
the noise added by the amplifier can be the crucial ele- 
ment in determining the feasibility of a correlation func- 
tion experiment. It becomes even more important for 
higher order correlations, where the number of samples 
depends on the noise power raised to some larger expo- 
nent. 



Appendix B: Coherence functions for states with at 
most one photon 

In the experiments described here [T] , the cavity is pe- 
riodically prepared in the state a\0) + with a period 
tp such that Ktp ^ 1. This ensures that, to a very good 
approximation, the cavity returns to the vacuum state 
before the superposition is prepared again. 



The coherence functions can be calculated straightfor- 
wardly via their definitions in terms of the field correla- 
tions, while the correlations can be calculated by solving 
the Heisenberg equations of motion for the cavity field, 
and using the quantum regression theorem (331 Ell HZ] ■ 
This procedure can be greatly simplified by noting that, 
if t and t + t are in different preparation periods, then 

{a\t)a{t + T)) = {a\t)){a{t + T)), (Bl) 

and 

{aHt)a\t + T)a{t + T)a(t)) = 

{aHt)a{t)){aHt + T)a{t + T)), (B2) 

due to the assumption ntp ^ 1. 

In the case were t and t+r are between ktp and {k-\-\)tp 
for some integer fc, one finds that 

{d\t)a{t + t)) = (n(0))e-"(*-'='p-"/2\ (B3) 
{a^it)a^{t + T)a{t + T)a{t)) = (a^ataa)e-«(2*-2fctp+-), 

(B4) 

while if t and t + t are in different preparation periods 
starting at ktp and (fc 4- l)tp, one finds that 

{a\t)a{t + t)> = |(a(0))|2e-'^[*-'=*p-('-"p)/2l, (B5) 

{a\t)aHt + T)a{t + T)a{t)) = (7i(0))^e^''(2*-2'=*p+^-"p). 

(B6) 

After integration over t, the first-order coherence func- 
tion can be shown to be well approximated by 

G(i)(r) = -(n(0))e-«l^l/2 

K 

+ -I(a(0))|' Ve-"l"-"pl/^ (B7) 

1^0 

This can be interpreted as a series of time-shifted copies 
of e~''l'^l/2^ where the peak centered at r = has a height 
equal to (n(0)), while the peaks centered at non-zero mul- 
tiples of tp have a height equal to |(a(0))p. 

Under similar assumption, the second order correlation 
function can be shown to be well approximated by 

G(2)(t) ^ -(a^(0)a^(0)a(0)a(0))e-"l"l 

K 

+ -(n(0))'y e-«l"-"pl, (B8) 

such that the center peak has a height proportional to 
(d^(0)a^(0)a(0)a(0)) while the other peaks have heights 
proportional to (n(0))^. 

For the superpositions of vacuum and a single photon 
considered in [T1, we find that 

(n(0)) = WW (B9) 

|(a(0))p = |an/3p, (BIO) 

(a^(0)a^(0)a(0)a(0)) ==0, (Bll) 



(B12) 
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indicating that the center peak of G*^^^ is abscent, while 
the other peaks are non-zero, which is a signature of the 
purely quantum effect known as anti-bunching [71 |H] • 



Appendix C: Filtering 

The finite bandwidth of the detection chain can be 
modeled by considering the insertion of a bandpass filter 
in an ideal (infinite bandwidth) detection chain. In order 
to calculate the effect of filtering on correlation functions 
one can consider a general framework which describes 
what happens to multi-time, multi-channel correlations 
when measurement signals are filtered. Assume a system 
with n channels where each channel is filtered individu- 
ally. One can write the filtered outcome of each channel 
Sfii,i in terms of the input signal Si and the filter func- 
tion fi by using the relations for linear time-invariant 
systems [48] 

Sfii,iiti) — fi{ti) * Si{ti) 

n+oo (CI) 

fi{n)Si{ti - Ti)dTi. 

Each channel has a separate time variable ti to capture 
the case of multi-time correlations. This also clarifies 
with respect to which variable the convolution is done. 
The goal is now to express the filtered coherence function 

Gfii(ii, . . . ,tn) = (<S'fil,l(tl) 5'fii^2(^2) • • • 'S'fiLn(tn)), (C2) 

in terms of the unfiltcred coherence function 



G(tl,...,t„) = (5i(ii)52(i2)---5„(^„)). 



(C3) 



This can be done straightforwardly by substituting 
Eq. (ICll) into Eq. 



Gm(ti,...,i„) = i'Wh{U)*S,{ti) 



(C4) 



\i=l 



Realizing that all convolutions are related to different 
time variables one can rearrange this expression as 

Gfii(ii, . . . , t„) = 

flih) * /2(i2) * • • • fnitn) * G(ii, . . . ,t„). (C5) 



one can write 

Gmih, . . . ,tn) = F{h, . . . ,tn) * G{ti, . . . ,tn). (C8) 

In frequency domain, the same fact can be expressed by 
using the multi-dimensional Fourier transform instead, so 
that one may simply write 

Gmi^i, ■■■,^n) = G(a;i, . . . , w„)i^(wi, . . . , w„). (C9) 

1. Two-point correlation functions 

Using the spectral representation of some first-order 
coherence function G{ti,t2) and the global filter function 
F{ti,t2), one can write Gfnir) as 

Gfil(T)= / dtdt^idw2e^^"^+"^^*+'"^"i^(wi,w2)GK,c^2), 

(CIO) 

= — / duje"^''F{^uj,uj)G{-uj,uj). (CIl) 

Considering the time representation of this expression, it 
is clear that the correlation function will be distorted by a 
convolution with the effective two-point correlation func- 
tion /cff(T) = J^{F{—uj,uj)/2tt}. Due to the linearity of 
the filters, one finds that Dirac S in the noise correlations 
are replaced by /off, so that, for example 



(CI2) 



where h^c\t) have been introduced in Eq. (25) and (•)fj[ 
indicates that the average is taken over filtered outputs. 
This illustrates why the values for the different second 
order coherence functions remain finite. Moreover, the 
other time-integrated two-point correlations are replaced 
by the convolution of the two-point correlation function 
with the /off. 

Note that since a linear time independent filter is used, 
the relative heights of the peaks remain unchanged - only 
their shape gets distorted and scaled. This is illustrated 
in Fig. [6] and Fig.[7| 



The integral form clarifies this expression 



2. Four-point correlation functions 



Gfii(ii, . . . , t„) = 

/+00 /•+00 
dTi--- dTn /i(ii-ri) • • • /„(i„-r„)G(Ti, . . . , r„). 
-oo — oo 

(C6) 

This expression can be seen as a generalized convolution 
with respect to more than one time variable. Introducing 
the global filter function 

F(ti,...,t„) = /i(ii)/2(t2)---/„(i„), (C7) 



Considering the second-order coherence function with 
filtered signals, one obtains 

Cm (^) = J dtdoji--- da;4e*("i+"=+'^^+'^'')*e*('^^+"'')^ 

F(a;i,W2,a;3,a;4)G'^^'(a;i,a;2,a;3, W4). (CIS) 

The different types of correlations discussed simply de- 
termine the labeling of the variables. Applying a change 
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FIG. 6: (Color online) Distortion of the 11 pulses in G'^'(r) 
due to Gaussian filters of different bandwidths. The solid blue 
line is the unfiltered function, and the others are filtered band- 
width decreasing progressively: 31/tp for the dashed purple 
line, 14/fp for the dot-dashed yellow line, and 10/tp for the 
dotted green line. 




FIG. 7: (Color onhne) The pulse train for unfiltered G^^'{t) 
(top) and Ggj' (r) after Gaussian filtering (bottom) under the 
periodic preparation of \/V3|0) + \/V^|l)- The Gaussian filter 
used here has a bandwidth of 31/ip. 



of variables and integrating over time, one obtains 



1 p+oc 

(t) = — / dn2dn3dQie-'"^^ 



F 



'^2 + ^2+ ^3 ^2 - -^2 - ^^4 



2 ' 2 ' 2 ' 2 

-^2 + ^^4 + ^2 - ^^3 -f^2 - ^^4 



(CM) 



Note that this leads different behavior, in the sense that 
the correlation function is not simply convolved with an 
effective impulse response. 



Appendix D: Commutation relations of two-sided 
cavity outputs 

Taking both cavity mirrors to be leaky, one finds an 
additional boundary condition in the input-output de- 
scription of the cavity O [3TH33] 



(Dl) 



where the Cin,out mode are now the modes coupling to 
the second leaky mirror. The equation of motion Eq. ([9]) 
for a in the rotating frame then becomes 



; Kft + Kc ^ , I f , I 

a = a + WKbOin + VKcCu 



(D2) 



From Eqs. ^ and ( |dT| ), one finds 



- V^b[m,cl{t')] - y^c[bin{t),a\t% (D3) 

where the input field operators were taken to commute. 
Integrating the solution for the equations of motion of 
the modes b{uj,t) of the left transmission line and the 
modes c(w, t) of the right transmission line, and using 
the definition of the input fields one obtains 



+00 

dio b{(jj, t), 

(D4) 
duj c(w, t). 



Cin(t) 



a(t) 



1 



From causality and the boundary conditions above, one 
finds [53] 

[a(i), 6in(t')] = 0, W). Ci„(t')] =0 for t' > t 



(D5) 

[a(i), 6out(t')] = 0, [a(i), Cout(t')] =0 for t' < t 

(D6) 
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Combining these commutation relations with the input and therefore 
field definitions, one finally finds 

m,cl{t')] = ^Mt^t'Mt),aHt% ^^^^ 

[binit),a\t')] = y/^,u{t' - t)[a{t),d\t% [^outW>cLt(i')] = 0, (D9) 

where 

r 1 t>o, 

u{t) = < I t = 0, (D8) as claimed. 

[ i < 0, 
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